import control  #pip install slycot  #pip install control
import control.matlab

#unit s
Ts = 0.001 
#sys = control.TransferFunction(5.235e005,[1,87.35,1.047e004,0])
sys = control.TransferFunction(400,[1,50,0])
dsys = control.matlab.c2d(sys,Ts,method='zoh')
(num,den) = control.tfdata(dsys)
print("numerator = %s" % ("".join(str(lst) for lst in num)))
print("denominator = %s" % ("".join(str(lst) for lst in den)))
#numerator = [array([0.00019671, 0.00019346])]
#denominator = [array([ 1.        , -1.95122942,  0.95122942])]
